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(54) Abstract Title 

Procosaing of afgnals inecdont on an array 

(57) A method of signal processing comprises sampling a 
ofements of an array and transforming these signals betwec 
time delay to each of the transformed signals to align the 
applied to each of these signals for beamsteerlng and the 
inverse transformed between time and frequency domain tci 
signal. The method is carried out using a signal processing 
(4) for transforming signals Incident on the array between 
a time delay to each of the transformed signals to align the 
applying a phase shift to the aligned signals; a summer (7) 
and an inverse transform (6) for transforming the summed 
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Incident on a plurality of 
main, then applying a 
parate phase shift is 
ignals are summed, then 
and correlated output 
an array (1), a transform 
; a delay (5) for applying 
phase shifter (6) for 
fned transformed signals, 
frequency domain. 
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METHOD AND APPARATUS 



n 



This invention relates to a method of $h 
wide band signals. 

Particular applications of signal procesislng include 3Ct\jb rfttmote sensing, synthetic 
aperture processing and spread spectrum conmuntcatton in t^ e f i Ids of sonar, radar, and 
space communications. There are various known methods of Lir ring out the processing 



steps. Signals incident perpendicular to a stra 
Signals which are incident at an angle will havle a different ph^ 
array, but the array can be steered electronics liy to compansat a 
beamforming, an electronic time delay is applijsd to the signal ff^n* 
them to be summed coherently. In frequency 
antiy are transformed to the spatial frequency 



domain in which 

wave con^esponds to the spatial frequency. This transformatior 



applying a Fast Fourier Transform (FFT) to the 



Frequency domain beamforming is prererable because 



rate, the Nyquist frequency, two samples per cycle, compared y vtti' I the multiple samples 



needed for time domain beamfonming. However, a problem of Usi 
they only work for narrowband signals. For vAiie band signals, me 
depends on the frequency. There is a requirement for the syst m i 



process wide bandw'dth pulses to achieve short range resotutic f i i lid low reverberation 



FOR SIGNAL P ^O 



gnal processing 



set of array sigial ; 



levels. In sonar systems the signal reflected from the surface, 
temied reverberation; it masks low target strength echoes; this 
area insoniTied and thus to the range resolution. A short pulse 
however, only has a small energy and can be affected by noiS€ 



amplifiers in the equipment or from rain. Comrfionly, a long, vm*< ht wnd pulse is transmitted 



and the received pulse is compressed by correlating it v^h the 



very computationally intensive and the process ing required ina sas i^s with the bandwMth. 



In accordance with a first aspect of the 



present inventior t 



processing comprises sampling a time sequence of signals ind i ^r 1 on a plurality of 



elements of an array; transfonning these sign^ds between time 
applying a time delay to each of the transforms id signals to aligili 



applying separate phase shifts to each of thesis signals for beai r st i ering; summing the time 



aligned transformed signals; and inverse trans 



and frequency domain; whereby a beamformed and correlated 



ght line an^y m^jf t i summed coherently. 

|t each element of the 

f<! 



this, so in time domain 
each array to enable 
jomain beamfor||i|iii|||i, the signals from each 
\ ha skew in the incident sine 
is Itisually achieved by 



;i ^a : ed and water column is 
^ro| : lem is related to the 
\ good range rasolution), 



^Vll 



»OL rces, such as from pre- 



lar emitted pulse. This is 



ni^fuires a lower sampling 



; articular for receiving 
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fg spatial FFT's is that 
angle of electronic steer 
described to be able to 



Tiethod of signal 



nc frequency domain; 
Ui I signals for correlation; 



orming the sumr \bc ' output between time 



ut signal is produced. 



able to produce 



l)iiamfonneci and 
; gn als. The sampled signals 



up table. 



In accordance with a second aspect of t ie present invenfi 



The method of the present invention U 
correlated output signal by summing time aligned transfonned 
are derived from a transmitted signal which is 
sections, each sub-secton comprising an indi pendent band o1 
bandwidth less than the bandwidth of the signal and wherein e ich sub-section has a 
duration greater than or equal to the inverse of its bandwidth 

Typically, the bandwidth b of each sub section is less tfi s n it >t 
bandwidth B divided by the square root of the 
transmitted pulse length, T, i.e. b <= B / (BTf^ 

Typically, the duration t of each sub-se :tion is greater tt i m 
duration T divided by the square root of the pn>duct of the sign; i 
transmitted pulse length, T, i.e. t T / (BT) 

Preferably, the signals are sampled at i\ sampling rate o 

Preferably, the time delay applied to each signal is an in 
of each transfonned signal. 

Preferably, the signals comphse quantised chirp signals 

Preferably, the component frequencies of the quantised 
harmonically related and have equal durations. 

Preferably, alternate sinusoidal components of the aligns <l $ gnais are converted into 
cosinusoidal components by applying a suitabli^ phase shift, wh i 
peak is generated when the final output is generated by Fourier 



Preferably, the phase shifts for beamstcering are precalc i fat »d and stored in a look 



constrained to a mp^ise discrete sub* 

re<|uency components of 



equal to the total 
product of the sil|)al bandwidth and the 



i;i6r equal to the total 
bandwidth and the 



twt:6 the signal bandwidth. 
i r multiple of the duratiori 



[hir> signals are 



refcy a single congelation 
syrihesis. 



device comprises an array; transfomr> means fo ' transforming sigr ia|i ; Incident on the array 



tirr 



e delay to each of the 
rpeans for applying a 

the time aligned 
t >e summed output 



ng 



between time and frequency domain; delay me^ins for applying 
transfonned signals to align the signals for com Nation; phase sh 
phase shift to the aligned signals; and summing 
transformed signals, and inverse transform me£ ns for transformi ig 
between time and frequency domain; whereby n beamformed ar ( I c( stated output signal 
is produced 

An example of a method of signal proM^^sIng in accordar (a 
invention will now be described with reference t<» the accompany r ig 

Fig. 1 is a block diagram for apparatus f(f frequency dom a 
according to the present inventton; 



a signal processing 



Mth the present 
itirawings in which:- 
ini! signal processing 



Fig. 2 is a S-dimensional plot of a signal processed in a|fojf|iance with the present 
invention; and. 

Fig. 3 is a block diagram for apparatus 
to the invention. 

The apparatus described in this example is for signal prlblc^sing from a sonar array, 



although it is equally applicable to communica tk>ns and radar. 



for time domain 



operates in the frequency domain and shows 



I sonar array 1 vi rirc i receives signals derived 



from a transmitted waveform reflected from ot jects onto each if diiient of the array. The 



systems describ^q ij 

6 



preferred transmitted vraveform for use in the 
waveform. K is a long pulse made up of a nurjiber of sections A 
equal to one cycle of a fundamental frequency \ f. Each sectior 
wavefomn of different frequency, the frequenc 
be in ascending order. For a frequency v/hich 
frequency f then h*f cycles of the h th harmon 
smooth transitions from the end of one sectiori 

Although, the signal in each section is 
signal starts and stops after a particular time, ^. The bandwidtlli 

The total pulse 



esfonning a hanK)|i 
is an integer mut p|i 
c fit exacUy into o icj 
to the start of th 
"nominally of the lifarjii 

oi 



gi A\ processing according 



Ih. 



^n*s. i.e. n^ With a 
fn), a Fourier transform 



total bandwidth B of n adjacent sections is n/s 
n*s, so the BT product which is a measure of processing gain i 
quantised chirp waveform used for correlation (a matching wavls 
of size n separates the frequency from each section of the mat Jiipg waveform into a 
different independent frequency channel feed ng the delays 5 1 1 Ft 9. 1 . 

The signals from the array are basebaiided 2 to reduce I sampling rate needed, 
and sampled into a buffer 3. This operation s hifts the f requen ifh ^mponents of the 
received signals from the elements of the array, so that the hig lesi frequency in the carrier 
frequency chirp maps to the fundamental of \i e basebanded 1 rp , and lower harmonic : 
frequendes in the carrier chirp map to subsequent higher harmcinj;s in the basebanded 
chirp. Then the basebanded time samples of the signals from element are 



system of Fig. 1 



a quantised chirp 
ual duration, precisely 
c4ttains a sinusoidal 
!c series which need not 
h of the fundamental 
section. This allows 
rj|(xt, 

onic frequency h*f, the 
he signal is 1/s and the 
9r||tii, T of n sections is 



transformed as a block by a FFT 4 to generate complex numb^m 
each frequency component. Separate time delays 5 are applii 
coefficients from the FFT 4 to align the sinusoidal waveforms f^ 
synthesis of a the compressed waveform. A phase adjustmen 
from a look up table 7 to the frequency coefficients to rotate th 
adjustment coefficient is dependent upon whiiph element of the 
signal frequency and which beam, of plurality 
an area of interest, is to be formed. The phas 



of beams which c^rrji 



fi} quency coefficients for 
til the frequency 
r sl^bsequent Fourier 

Efficient 6 is applied 
r jnase. The phase 
irtl^y the signal is from, the 
an image representing 



e shifts are pre-c s lci . lated to correspond to 



4\ 



e 



ichieved in eithc r 



the delay needed for a simple time domisiln be imfomier The ctimf^lex 
coefficients are summed 8 over all the elements and input to a 
of an Inverse FFT 9 for the particular beam. Tpe output 10 of i 
time series of the data signals. 

Implementation of the method may be 
steps. In either case a large number of small i^FTs are parforr||i|ed 
shift registers or ring buffers to implement the iltgnment delays 
Accessing a large number of data stores can t e time oonsumir|( 
a separate pointer to access each, but this car 
single rfng buffer array with a set of interlinked 



processed signal 10 is output to a display through a graphical ii i 

ocks (numbered 



be overcome b; ^ 
pointers to conti ^ 



nd 



Fig. 2 shows a 3-D plot of successive b 
processed output data signals. The time doms in output is plottiyd 
cof>erently processed output on a 16-sample a|ds (SI to S16), f 
blocks at successive instants of the array data 
in the orthogonal axis (numbered 1 to 65). Raijige sidelobes (1 
the beginning and end of each processing block, particulariy for 
waveform. These will manifest themselves on target echoes 
spaced at a sut>-pulse apart, and as lower amplitude echoes wHI6sji 
with the random offset of the position of the start of the process|r|)gi 
main target. 

It can be seen from Fig. 2 that the amplitude of the sideljbbo^ 
middle 14 of the processing block, but become^ significant at th||s 
the middle sectk)n of the processing blocks (i.e 
filling in the signal between adjacent processing blocks with a si 
processing blocks, the amplitude of these ranges sidelobes can 
level. Equally sized interleaved blocks give adc quate performai 
could improve the processing efficiency or the sidelobe suppresj 

Fourfer synthesis consists of summing s nusoid or cosin 
a given waveform. If cosines are used, two pes ks are generate 
processing btock and a second at the end. Pha 
coeffidents can be used to align the peaks to pioduce a matchejk 
response in the time domain. If sinusoids are u^ed, the summ4l< 
peak straddling the mid-point of the processing 
cosines by phase rotating (multiplying the comp 



>k>ck. The sinuspp 
ex frequency 



frequency 
appropriate frequency input 
f 'FT is a beamformed 



oir 



hsjrdware or software 
and the data stored in 
the correlation step, 
bti^cause of having to load 
irri >lementing all delays In a 
tf^ access. The 
er| ace. 

1 to 89) of coherently 
for successive blocks of 
successive processing 
) .! sampling rate are plotted 

13) appear towards 
1^ cosinusoidal source 
awlitional sidelot)es 
amplitude fluctuates 
[..block relative to the 



is quite low near the 
edges 15. By using just 
sigpial at the edges), and 
\xa id set of interteaved 
kfpi to an acceptable 
although other ratios 



on 

id$ 



on. 

I^oii I waveforms to produce: 
-cneatthe start of the 
;^m^te (complex) frequency 
ered output of a single 
produces a double 
can be converted into 
co|ii(ip4nents by y ), and then 



applying the phase reversal technique to gen 
filtered waveform. 

The complex coefficients needed to produce the phasd 



hamnonic and beam combination are preferably stored as a lo( \ (u|; 



changes needed to align the peaks for sine o^ cosine matchint i 



incorporated into the coefficients stored in the 
Correct alignment of frequency bands 



section as received from the sea bed, for corr diation by Fourie ' sjf ithesis can be checked 



by displaying or printing the complex frequency components c i 
processing. The simplest test is to apply an impulsive echo wt t::H 
chirp back to the receiver. After n blocks hav^ been processejl (ri 



frequency components of a particular beam should have their na): mum amplitude 



simultaneously. Numeric display methods are 

multhbeam data, since the method of signal ph)ces5ing operat|^ wide band signals 
consisting of many samples, so a graphical irrerface is used 

In a second example operating in the t me domain shov|f^ '4 
each element of the array is transformed by d claying it by one 
and subtracting in a subtracter 18. passing the output through 
separating the frequertcy components into separate time doma 
variable digital filters, otherwise known as res9nators 20, to ini^^td 
terminate it after exactly one sub^ulse. 

This operation provides time samples 6f the signal sep^rjati^d into its component 
frequencies, and requires sampling at least twjice the Nyquist r^e 
bandwidth of the data signal. However, even 
vernier time delay for multi-element beamfomntng. For exampi j 
beamformer for a 32-element, 32-beam system with a bandwid J[t < 
order of 5 MHz sampling rate, whereas twice riyquist is only 25 



»rate a single ce 



I f^i I peak to the matched 



table, 
comesponding tjb 



not appropriate 



i:ri lecking time domain or 



Fig. 3, data 16 from 
|s|ub|pulse in delay line 17 
tapped delay 19 and 
ifi^gnals using state 

a resonant response and 



both in-phase (Year) and quadrature (Imaginary") time domain 



of two integrators which comprise the resonate ^r stage 20. The » M f e used to modify the 
phase of the resonator output by a small angUf needed to align \^ 
the incoming wavefront at a desired steer ang e. The alignmer 
components • a coarse delay 22 provided by an offset in a ring 
component 23 produced by phase rotation usipg the complex c(|::|m|^nents of the time 
domain signals. 



It is however possible to provide vernier adjustment sine ( tiie resonator provides 



sh fts for each element. 



wiivefomns maybe 



friiquendes in each sut>* 



table 10. The phase 



Clip ated at each block of the 
reflects the transmitted 
-sqrt BT), all the 



e. four times the 
00 tow to provide a 
create a synchronous 
64 kHz requires of the 

kHz. 
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sari^ples 21 from the outputs 



signals from elements to 
thin consists of two 
bufler, and a fine or vernier 



It might be thought that the individual n^sonator outputs 



frequently than twice Nyqulst of the wide band 
(1/sqrt(BT)) of the total bandwidth of the signa 



signal, since the 
. However, this 



accurately because of the out of band comporants which pass 



resonator components. The time signal ateacTi 



frequency, and in the frequency domain this maps to a sync fur (^Ofi rather than a single 
frequency. Measurement of the magnitude of |he complex vec^r spows that its amplitude 
is modulated within the sub pulse period when 



the time domain 



fed with the matching waveform, and samples faken at Nyqulst 

subject to simila 



c O0T elator resonators are 
i>rtie sub-band (i.e. 

1/sqrt(BT) of the overall Nyqulst rate) would be subject to simila flMi:tuat'ons because the 
sampling instant for the lower sampling rate ocfurs randomly re|4tiv^ to the position of the 
peak in the echo retum. 

The phase shifted frequency coefTicientk aligned for confell 



s tew ed filtering and 



across all frequencies, and the resulting compcnents of the surr rne^l beams are summed 
25 for each element, then input to an inverse FfT 26 to genera^jd ij series of time domain 
outputs 27. 

The method of the present invention carji be applied to n^|g 
beamforming processing of wide band signals, i 

including synthetic aperture processing, and spfead spectrum a|i|nfTjunication, in fields 
including sonar, radar, airt>ome sonar, and space communicatio 
Additionally, since the intenmediate frequency-separated 



signals are directly controllable within the process, many additioi i^il { recesses can be 



o tr e capability of the 
< )f I ndivkiual frequency 
{:3n >e used to increase the 
/|e» >lution of the sonar to 
he spectral spread of 



incorporated into a beamforming/correlation processor to increa< 
overall package. Examples are separate detection (rectification} 
bands prior to signal combination In active sonap systems. This 
revert)eration to shadow ratio or alternatively to 
the expected target size.. Another application \i in measuremen||of 
the retum signal which provides a measure distipguishing point r||n^ors from random 
reflectors. 

Processing efficient using this method 
from time domain or FFT beamforming techniques, thus allowing 
with greater processing gain against ncMse and r sverberation. In ^ 
separately steered, to provide high resolution wi idows within a Ic v^r 
or simply an electronically steerable window. Bcjam focusing, anfi 
be achieved allowing operation in the near fiekJ 



qou d be sampled less 

ei ich comprise only 
^bpioach fails to operate 
itti^nuated, through the 



resonator is a t i\e- imited pulse of a single 



15. 



End 



of a conventiona 



on are summed 24 



remote sensing. 



spatiafty-separated 



i:ar\ be greater tfjsn |i;onventiona!ly available 



US| 



of wider BT products 
victual beams may be 

resolution backdrop, 
diyjnamic focusing, may 
arr ay or of a synthetic 



array produced by the motion of a real array. 
beamforming algorithms and height finding algorithms. 



Other appiicatioi ii> s e high resolution 
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signals incident on a plurality of elements of 
time and frequency domain; applying a time 



1 . A method of signal processing; the vr ethod comprisim i sa npiing a time sequence of 



an array; transfo inir g these signals between 



align the signals for correlation; applying separate phase shif s to 



beamsteerfng; summbig the time aligned transfomried signals 
summed output between time and frequency domain; wheredj^ 
correlated output signal is produced 

2. A method according to claim 1, wherdin the signals an 
of twice the signal bandwidth 



clain 



3. A method according to claim 1 or 
signal is an integer multiple of the duration of 

4. A method according to any preceding 
quantised chirp signals. 



2, wherein the 
each transfonmetk 



daim, wherein t\ 



5. A method according to claim 4, wherelp 
chirp signals are harmonically related and havie 



6. A method according to claim 1 , wherei i 
aligned signals are converted into cosinusoidsl 
shift, whereby a single correlation peak is generated 



7. A method according to any preceding 
precalculated and stored in a look up table. 



ar;iy; 



8. A signal processing device; the device 
transforming signals incident on the an^y betv^een 
means for applying a time delay to each of the 
correlation; phase shift means for applying a 
summing means for summing the time aligned Itransformed sigrf^ls 
means for transforming the summed output between time and 
a beamformed and correlated output signal is (produced. 



elay to each of 



lie lansformed signals to 



t^m^ delay applied to each 
signal. 



component fre i 
equal compon » 



alternate sinusjblda 
components b> 



each of these signals for 



mi inverse transforming the 
al>eamfonned and 



sanpled at a sampling rate 



s sgnals comprise 



claim* wherein thi \ 



:x)mprising an 
time and 
transformed si 
phase shift to the 



jemdes of the quantised 
itlsngths. 



components of the 
sipf^lying a suitable phase 



phase shifts are 



transform means for 
fre||:}ue|icy domain; delay 

to align the signals for 
a^igrjied signals; and 

and inverse transform 
f|i^qiiency domain; wtiereby 
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